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Part I 

SLOW GROWTH OF CRACKS IN A RATE SENSITIVE TRESCA SOLID 

by 

Michael P. Wnuk 
(Abstract) 

The paper proposes an extension of the classical theory of fracture 
to viscoelastic and elastic-plastic materials in which the plasticity 
effects are confined to a narrow band encompassing the crack front. 

It is suggested that the Grif fith-Irwin criterion of fracture, which 
requires that the energy release rate computed for a given boundary value 
problem equals the critical threshold, ought to be replaced by a differential 
equation governing the slow growth of a crack prior to the onset of rapid 
propagation. A new term which enters the equation of motion in the dissipative 
media is proportional to the energy lost within the end sections of the 
crack, and thus it reflects the extent of inelastic behavior of a solid. 

A concept of "apparent" surface energy is introduced to account for 
the geometry dependent and the rate dependent phenomena which influence 
toughness of an inelastic solid. 

Three hypotheses regarding the condition for frac.ture in the sub critical 
range of load are compared. These are: (a) constant fracture energy 
(Cherepanov), (b) constant opening displacement at instability (Morozov) 
and (c) final stretch criterion (Wnuk). 
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Slew Growth of Cracks in a Rate Sensitive Tresca Solid, 


by 

Michael P. Wnuk 

Theoretical analysis shows that the amount of slow growth occurring 
in a plane stress tensile specimen, subjected to a subcritical stress intensity 
level, is affected by the following parameters: 

(1) ductility and rheological sensitivity of the material, 

(2) rate of loading, 

(3) initial crack size, 

(4) geometrical configuration of the test. 

Some of these factors were included in the earlier serai-empirical approaches 
of Krafft et al. [3], Brown and Srawley [1], and in the only available exact 
treatment of slow growth under tearing mode of deformation due to McClintock 
[4], McClintock and Irwin [5] and Rice [6], here Rice’s idea of a universal 
R-curve is re-examined. It turns out that the curve is universal, that is 
independent of geometry and the initial crack size, only in the limiting case 
of a rate-insensitive elastic-plastic solid. When time effects are accounted 
for in the constitutive equations of the matrix which contains the crack, the 
"universality" no longer holds. 

The governing equation which describes the quasi-static extension of the 
crack under confined yielding condition is an integro-dif ferential equation 
derived in. {8]: 


R - 


- 4) + 4 Jr + f u 

+ f r “ 2 /<y? <SHP(S> - a) t ili : • 


(1) 
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Here R denotes the length of plastic zone ahead of the crack front, R q 
is its value at the growth initiation, A is the dimension of the process 
zone, ¥ is the normalized creep compliance. Both R and Y are functions of 
time t, or equivalently, functions of the current crack length l, which is 
treated as a time- like parameter. Other relations are 


£ = A/R(t - fit) , fit = A./& , i - dZ/dt 
p (s ) = 1 + (e - s) dR/d£ , 6V = Y(fit) - Y(o) 


The criterion for crack opening, which was employed in order to derive 
eq. (1), is that of "final stretch", cf. [8]. In contrast to the COD 
criterion the final stretch condition is path-dependent and thus it appears 
to withstand Rice's criticism [7] of earlier work on this subject by Cherepanov 
[2]. Postulating this criterion in [8] we required that the amount of stretching 
vhich occurs within the process zone during the time interval just prior to 
fracture ia a material constant. Such an approach assumes nothing about 
the current tip displacement and the length of the associated plastic zone. 

In fact these two entities turn out to be functions of time and the loading 
history. 

It should be noted that the final stretch criterion coincides with 
McClintock's criterion of critical strain attained over the Neuber domain. 


One may also addjin the limit case of steady-state propagation, both the 
COD and the final stretch criteria converge. The essential difference 
between them becomes obvious, though, within the subcritical range of 
applied stress intensity. 


To ijiake the problem mathematically tractable we assume further that 
the length of the process zone is small vs. the plastic zone size, i.e. 
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A/R <<1. Then eq. (1) reduces to 

A 4R -1 dR -1 3R 

(2) j •'M'S”) + A • [1 + CROR/3Q) ] - AC(9R/9Q) = R Q 


Here (dR/dft- 9R/9£)CA (9R/9Q) ^ has been substituted for ^^(“BA/i) ; B denotes 
the slope of creep compliance at time zero, Q is the loading parameter and 
C = B/Q. To illustrate applications of eq. (2) we shall integrate it 
numerically for the case of a crack contained in an infinite plate under 
tensions; then R =(l/2)Q^(£)£ , Q = Tr<a/2Y. Since the ratios R/A and It, /A 
are very large numbers, it is convenient to introduce the logarithmic function 


and to cast the eq. (2) into the following form 
dY _ A - (1/2)Y+ ^ exp (3/2) (Y - X) 


i + A 


‘JY exp (1/2) (Y -X) 


where 


X = log (ft /A) , A = (R/A) - log2 

(4) ° 

Y = log(R/A), Y = Y (X) 

Equation (3) has been integrated numerically for a certain set of initial 
conditions, and the results are shown in Fig, 1. The figure illustrates the 
effect of rate sensitivity of the material and the rate of loading on the 
shape of the R- curve • It is seen that not only the slope of the curve is 
affected, but also pronounced changes in location of the ultimate instability 
point are observed. Examples of integration of the equation of motion when 
the visco-elas tic dissipation is dominant, are shown in Fig, 2 and Fig, 3, 
Such a case of a n creeping crack" is considered in more detail in [9]* 

For engineering applications it is convenient to re-write the governing 
equation (2) in terms of the ratios R/R^ and fc/R*, , where denotes the 

2 o 

steady-state limit of the plastic zone size, say - ttK^/SY . The symbol 



4 


K ot denotes the maximum plane stress fracture toughness which would be attained 
in an ideal case, when the conditions of the test are such that the prior to 
failure slow growth is fully developed. Of course, the actual fracture 
toughness K^, i.e. the value of K at which the rapid motion begins, is bracketed 
by the initiation toughness and the maximum steady-state toughness K^, 
Normalizing the plastic zone size and the crack length as follows 

(5) R/R, = p , 1/ R^ = C , R ro = <;/8Y 2 


changes eq. (2) into 


(6a) ■§ - 
or 

(6b) ija = 


(l/2)g.n(i/p) + Cp 2 /;: JTEp 

i + c p/yw 

fcn(2/0 2 c) - Q 2 
2CQ(1 + CQ/2) 


p =(1/^ 2 C 


for a crack in an infinite plate, and 

(J ) _ (l/2Un q/p) + Cp 2 {l + (pe/co)tan(pC/Co)}/gy2pgsec(pc/e^)' 

^ 1 * Cp/ /2p tsec(p^/ CoJ 

or 

p = (l/2)} 2 £sec(ps/co) 

dQ _ ^,nt2/Q 2 Csec(pc/co) ] -Q 2 [l+ (pe/co) tan( P C/co) Isecfpt /g qj 

2Qt(l+ ^Q/2)sec(p^/to) 

for a crack traversing a panel of width 2K The initial crack length is given 
by = &o/Rco and P denotes the initial crack length to panel width ratio , 
p = Tr£ Q /2b. The locus of terminal instability follows readily from eqs . 
(6b) and (7b) if dQ/dl is set equal to zero. Then for a crack contained in 
an infinite plate one has 

(8) C f = 2/Q 2 • exp (Q*) 



while for a central crack in a finite width panel it is 

(9) in [2/Q 2 f C f sec(pc f /C 0 )l = Q 2 f [1 + (pC f / Z Q ) tan(pc £ /c 0 ) ]sec(pc f /c 0 ) 

Interestingly the rate sensitivity C does not enter explicitly in the above 
relations. It is present here, though, in an implicit way, since both the 
critical load and the critical crack size are pronouncely affected 
by the rate sensitivity. This can be seen only after the integration of 
equations (6) and (7) is completed. See Fig. 4. 

The effect of finite width on the amount of slow growth which takes 
place prior to failure is illustrated in Fig. 4, where Q vs. c curves are 
shown for both infinite and finite plates at certain levels of loading rate 
Q. Two different trends are observed: (1) the slow growth is enhanced when 
the panel width stays - constant while the initial crack size is increased, 
and (2) the slow growth is diminished when initial crack size Is kept constant 
but the panel width increases, see Figs. 5 and 6 . 

In general the amount of slow growth before the final instability sets 
in turns out to be a function of (1) ductility, (2) rate-sensitivity, (3) 
rate of loading (4) initial crack size, and (5) geometry of the test. Although 
no closed form solutions are available at this time, the influence of the fore- 
going factors has been investigated numerically, and the results are gathered 
in Figs. 5 through 8. The graphs were derived from a number of integrations 
of the governing eq. (2) performed on an IBM 360. 

It should be noted that for the case of fast loading, or equivalently, 
for a rate insensitive solid (C-H3) , the equation (2) which describes the 


R-curve degenerates to 
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This has a closed form solution 

(11) £ -Z 0 = 2{ei[£nj> 0 ] - ei[ £ n £ ] j 

x t 

ei(x) ~ f dt 

•— OQ 

It is seen that the last equation supplies a universal relationship* since 
neither geometry nor the initial crack size are represented (the initial crack 

length enters only in form of the difference C - C ). Thus the shape of the 

o 

resulting R vs . £ curve will remain unaltered by these factors. On the other 
hand* the location of the terminal instability point will depend on a specific 
geometry of the test. It is so, because the instability state is determined 
by the point of tangential contact between the R-curve and a member of the 
family of curves originating at (R=0,£=0) and representing the variation of 
plastic zone size with crack length at certain fixed values of loading parameter . 
Slope and the shape of these curves will, of course, depend on the geometry 
of the problem, and so will the position of the instability point. 
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APPENDIX 

A. Fracture Criterion and Slow Growth of a Crack in an Elastic-plastic Solid. 

A local criterion of fracture is postulated by proposing that the sum 
of work done at a fixed material point, which undergoes a deformation process 
while it traverses the Neuber section of the plastic zone, is a material 
property. Thus the crack will move onward if 
t 

(Al) / S[x ,t]u[x . ,t] dx = <f> A 

t-<5 1 F * 

Here, S[x^,x] is the restraining stress at the control point P and time x, 
while u [Xp,x] denotes the time-rate of the displacement at the point P 
and time x, and perpendicular to the crack plane. Time increment 6t equals 
the time used by the crack front to pass through the Neuber zone (of 
characteristic length A) and thus it is related to the rate of crack growth 
as follows 

(A2) <5 1 = A (i 

Such a relation was used earlier by Glennie and Willis [Al] to describe a 
piece-wise linear approximation of an accelerating crack. The material 
property <j>* has a dimension of energy and it can be related to the threshold 
fracture energy, i.e. the specific fracture energy at the onset of crack 
growth . 

We assume further that the restraining stress is constant (=y) over 
the Neuber domain and that the work <f>* can be expressed as a product of y 
and the initiation displacement u o , i.e. the tip displacement at which 
motion of the crack sets in. Under these assumptions criterion (Al) reduces 
to the "final stretch" condition which reads 
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(A3) u(x p ,t) - u(x p ,t - St) = u o , St = A/£ 

In other words the increment of deformation generated at the point P, just 
before an infinitesimal element located at this point collapses, should 
remain constant during the slow propagation stage. Thiscondition is not 
identical with the COD criterion since it allows for a variable tip 
displacement (note that there are no restrictions imposed on u(x p ,t), but 
only the increment, say Au p , as defined by the LHS of eq. (A3), is said to 
be a constant). Of course, when the motion attains the steady-state limit, 
i.e. when the length of the plastic zone R remains constant and the crack 
runs fast enough to justify the quasi-steady approximation of Glennie and 
Willis over the entire plastic zone, both criteria coincide. 

Interestingly, the final stretch criterion is identical with McClintock’s 
condition of critical strain attained over the Neuber domain, provided 
that one defines the strains within the plastic zone of a Dugdale crack as 
follows 


(A4) e(x ) * £ q + — {-[grad u(x^) |} , 0 * x * R 

o ^ 

P 

here is the yield strain, £f is the plastic component of the strain at 
fracture and u(x^) is the displacement within the Dugdale plastic zone. 

To show the equivalence of the final stretch criterion, as given by eq. (A3), 
and the McClintock critical strain concept, let us consider the case of a 
general in-plane loading mode (either mode I or II or both applied simultaneously). 
For such a case the displacement which results from the appropriate boundary 
value problem formulation is 


u(Xj> 





x i ^ i 

2 \l R <V -f R < x l - xi] J J 


(A5) 
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Note that the coordinate x^ is used here as a time-like parameter. Since the 
length of the plastic zone Is a function of time, it will therefore depend on 
in an a-priorL unknown fashion. To emphasize this point the symbol u(x^) will 
be replaced by u = u(x^,R(x^)). 

To apply the final stretch criterion we have to compute the displacements 
at the point P at time t and time t - 6t. These two instances correspond 
to = 0 and x^ - A, respectively. From (A5) we have 


(A6) 


u(t) - U( V 0) - ^ RCO) . H W« - (§ x ) - s > 

- u(Xj-A) - H { j[R(A)(R(A)-A)J - |^( {r(A) - 


Applying the fracture criterion (A3) we arrive at 


dR r a Jr+Ar - A] 

(A7) RfA-jJ - /[R(R-A)1 + # ^ = R o 

where all R f s are taken at = A (or at the time x = t-S t) . This is the 
governing equation of motion relevent to the slow propagation stage in an 
elastic-plastic solid, under the in-plane mode of loading (mode I or II*) . 
Equation (A7) describes the universal tf R-curve n , which was earlier discussed 
by Rice [6] for mode III. This curve is also equivalent to the ,f G~curve Tf 
introduced by Brown and Srawley in 1964, [1]. 

Note that the initial slope of the R-curve described by eq. (A7) 


(A3) 


A , l>,Rp _ n l + I - jJjsdA -11 

i li )J 2 l/v* +r[cn„/A) -u 


is remarkably similar to the initial slope of the R-curve predicted for mode III 
by Rice [6] 


(A9) 


, dR v modelll 
M{ ^ 


- 1 - An (R 0 /A) 


If we introduce a parameter a (which equals the ratio of the plastic strain at 


fracture to the yield strain e^/e o ) , then both (A8) and (A9) can be re-written 


as follows 
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(A10) 


( r~ : — . i JT+ o ‘ - Ja 

( i£) „ J a < a + U + 2 £njl + a r +-'{? 


(a) 

T a t ^ a 

a - in (1 + ct) (b) 

The graphs resulting from eqs, (a) and (b) above, and representing our solution 
and Rice's are shown in Fig* Al. They converge for both small and large 
values of the ductility parameter a. 

In fact,. within the range of practical applications, R will be much 
greater than A* Equation of motion (A7) reduces then to the simple form 

dR 


(All) A§.R 0 -mf) 


This form results also from McClintock's critical strain criterion of fracture, 
if one defines the strains as in (A4) and makes use of the condition 


(A12) 


e(x = A) = e_, where e = e + 

1 £ for 


To demonstrate this point let us combine (A4) and (A12) 


(A13) 


grad u(x.) = u /A 

1 . o 

x * A 

1 


and let us compute the gradient of the displacement within the Dugdale 
plastic zone. Since u = u(x^,R(x^) ) , we have 


(A14) 


grad u(x ) |= T + ~ 
a 1 1 9R dx x 


where 


3u 

(A15) 'gg = 


4Y 

irE 


1 £ n ~ J&Kxj) ~ O ) 

2 + J( R ^ X P + x i)J 


dn 

dR 


41 r R C x l) - Xt , 

* E I l " ifcp 3 
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Note that dR/dXj can be replaced by -dR/d&, since x - l + x^ - Gonst. 
Substituting (A15) and (A14) into the critical strain fracture criterion 
(A13) , and remembering that R/A » 1, we get the equation of motion identical 
with the equation (All) which resulted from the final stretch condition* 

Eq. (All), describing the R-curve for an elastic-plastic solid, is a 
first order linear differential equation and it can be integrated in a closed 
form. The solution expressed in terms of A-units is 

(A16) - * 0 - (-|)jexp (2 R o /A)]{ei [An (fjs.) - 2^] - ei [An - ^2] } 

2 2 i 

while in terms of R. ( = ttK /8Y ) units it is^ 

K c 

R R 

(A17) A-A = 2R. {ei [An *2] - ei [An — ] } 

O w R* K* 

ei (x) = / [e /t] dt 

—CO 

Indeed, we observe that the amount of slow growth predicted by eq. (A17) 
does not depend on either the geometrical configuration or on the initial 
crack length. It is, therefore, a universal relation analogous to Rice f s 
[6] result for mode III. 

+) To allow for the change from the A-units to the R^-units , or from the 
micro to the macro-units, and to satisfy the boundary condition at the 
steady-state limit R one has to incorporate another postulate regarding 

the "opening stretch", namely 

u o - ( !> I (4R * M > 

This transforms the equation (All) into the simple form 

dR 1 . /ri 
d£ ” 2 
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B. Deriviation of the governing equation of motion for a viscoelastic- 
plastic solid. 

Consider a viscoelastic matrix containing a Dugdale crack and described 
by the constitutive equations 

t 3e 

(Bl) s,, (t) = / G (t - x) — dx 

-00 1 BT 

/ 3 e (t) 

s (t) = / G (t - t) — dx 

-CO 1 ^ 


The displacement perpendicular to the crack plane ahead of the crack front 
is given by 

t 

(B2) u (x,t) - u° (x, t) + / ¥ (t - t) u° (x , t) dx 

y y t y 

o 

where u°(x,t) is the associated elastic solution to a given boundary value 
problem* The normalized creep function f ( t) can be readily related to the 
relaxation moduli G-^(t) and G^Ct) . The lower limit of the hereditary integral 
in (B2) denotes the time at which a given point (x,y=0) enters the plastic 
zone* Consider now a moving crack whose front approaches a stationary control 
point P located at x, see Fig. A2. For a general in-plane tensile loading 
mode equation (B2) reads 

U y (x,t) = J” [R(t)(R(t) + £(t) - x) ] + 


4 ^(t) - 


X 


£n 


(B3) 


( 


[R(t) + J~[R(t) + &(t) - x] 


■) 


+■ 


xf R (t) - J [R(t) + A( t ) - x] 

+ / » (t - T) t JirctXrCt) + 2(t) - *)] + * { [ R < t)+Mt) -1) > 

t ^ 2 \Jr(t) - J [R(t)+£(t)-x]/ 

o ^ 

where R denotes the current length of the plastic zone and is the distance 
measured from the crack tip- To apply the final stretch criterion of fracture 


dx 


we need to evaluate the difference 
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Au p = u(x ,t) - u(x p ,t - fit) = 


u°(x p ,t) - u°(x p ,t - St) + 

(B4) t , t-St . 

/ y (x) u°(x p ,t - t)dT - J T(t) u° (x p ,t -St - t) dt = 


Au° +/ [^(t-t)-f(t-6t-T)]u° (x ,t) dx + 

-6 r 


o 

t 

+ / V (t - x) u ( x 3 x ) dx 

t-St e 

The first term on the RHS of the above expression, Au°, is known from the 

P 

elastic-plastic analysis presented in the preceding section 


(B5) 



R + 



J"[R(R-A)] 



J r + Jtr - at \ 

fR - J[R - A] / 


The second term in the final form of (B4) can be shown to be proportional 
to the second time- derivative of the creep compliance Y and therefore it 
will be neglected as only the first derivative of Y is retained in the present 
analysis. The only restriction imposed is that the change in Y due to the 
shift of argument by 6 1 is small. This means that the function Y does not 
vary rapidly between time zero and fit, and that within this interval it can 
be approximated by just the first two terms of the McLaurin expansion 

(B6) Y(fi t) ~ Y(0) + Y(0) fit 


(of course ^ may vary within the interval (t Q ,t) but it is roughly constant 


inside each fit section). 
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The very last term of (B4) is further reduced as follows 


(B7) 


/ * (t - T)u°(x p ,x)dT = ~ ¥(0) } { J [R(T)(R(T) 7 l(t) + A(t))] 

t-6t t-6t 


2 



X TR(t) - Ut) + A(t) 1 
J[R(t) - Jl(t) + 


Note that the coordinate x^ is fixed and equal 9. (t). All functions appearing 
in the integrand of expression (B7) , although unknown, can be represented by 
the following Taylor expansions 

£(t) = £(t-6t) + Z • St 

«.(t) = £(t-St) + l • (x-t+6t) 

( B8 ) 

R(x) = R(t-5t) + ft • (x-t+6t) 

R(t) = R(t-6 1) + ft • «.t 


where both Z and R are considered constant within the 5t interval. Inserting 
the above expressions into (B7) produces 


<B9) 


f ^ (t - x) u° (Xp ,x) dx = 
t-5t r 


t * 

O’ t JZp-ctmpct) - 1 ;,V0 >i 


t -6 1 


_ i •( t ~ T ) An f JpCO ± JJrOO - £ ' (t— x) /R(t-<$ t) 1 \ \ A-r 
2R(t -St) W p C T ) - J[p(x) - £ (t-x)/R(t-6t)] /J 


Here p (t) denotes the ratio R (r)/R(t - St). Conibining eqs. (B9), (B5) 
and (B4) (subject to the final stretch criterion) gives the following equation 


of motion 
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(BIO) 


r dRA / J~R + J R-A 

R- J[R(R-4)] +A— 

iy „2 / ilr^s , s .7 s / Jp(s) + JIp(s) 
— R / ($>$(p(s) - s)J- - 


- s] 


p(s) - J [p (s) - s] 


-J } ds = R 


Note that 6 'J' = Y(6t) - ^(0) and the current time t has been replaced here 
by a dimensionless variable 


(Bll) 


l . 


s = 


(t - T) 


R(t - <S t) 


The upper limit of the integral in (BIO) is now defined as 


(B12) • ■ TuFTTt) 

while the function p(s) is given by 


(B13) ' p(s) - 1 + (e - s) dR/di 

Summarizing this section we may say that the function R=R(£) is determined 

by the non-linear integro-dif ferential equation (BIO) , subject to the initial 

condition R=R Q at & = ^ 0 * In order to make the problem tractable we shall 

reduce the equation (BIO) to a differential equation. This aim is achieved 

e 

by observing that for R/A » 1 the integrand of the integral / { } ds 

o 

can be expanded into a series 
(B14) { ) = 1 + 0(e) 

where 0(e) denotes terms containing e ( =A/R) and being on the higher order 
of magnitude. Retaining only the first term of (B14) , and neglecting A vs . R 
in the second and forth terms of (BIO), reduces the governing equation of 
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motion to the following form 

(B15) A df + 2 + 6f ' R = R o 

This is the desired form of the differential equation which allows for direct 
determination of the R function. Note that the LHS of eq. (B15) can also 
be written in a more compact form 

(B16) Aup = Au° + u° (Xp , t— 6 1) 6 'R 

Next, the increment SV is related to the R function. We have 
(B17) SY - ¥(6t) - 'R(O) * ’i'COjSt 

Now, denoting the material rate sensitivity f(o) by B, and recalling that 
it - A/£,we may further write 

(B18) &V = BA /£ - BA H - BA ^ (|£) “ 1 

But R depends on i in the following way 

(B19) R - R [A,Q(Jl)] 

where Q = Q(&) is a function describing the applied load (or dimensionless 
loading parameter). This generates 


(B20) 

hence 

(B21) 

and 

(B22) 


dR _ 9R 3R dQ 

dZ 3£ 3Q dJt 


is. _ /dR _ ML ,3Rv -1 
dZ ~ v d£ dZ } 



(45- 15>(2S) 

K dZ 


-1 


■Hit 

Let us denote the ratio of/^rate sensitivity B and the rate of loading Q by a 
letter C, and insert the result (B22) into the equation (B15) . This gives 
our governing equation of motion the following form 
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(B23) | + A[1 + CR (3R/3Q)" 1 ] J - 

- ACOR/3Q ) -1 R(3R/3Jl) - R q , C = B/Q 

Note that for C - 0, i.e. for a zero rate-sensitivity B or an infinite rate 
of loading Q, we recover the equation of motion valid, for an elastic- 
plastic solid, as shown in the preceding section by eq. (All). 

The other limit case follows from the equation 

(B24) Au° + u°(x , t - StJfiY = u 
p p o 

if one considers a purely visco-elastic solid. There plasticity effects 
are neglegible and one may think of A and R as being of the same order of 
magnitude (in other words the "plastic zone 1 ’ shrinks to just the "process 
zone 1 '). Of course A is still sufficiently small to justify the quasi- 
steady-state approximation, i.e. i * const, within the time interval of 
<St =A/f. Under these conditions the first term of eq. (B24) reduces 

(B25) Au£ - £ 7 * u°(t - fit) - I; u° (t - fit) = u° (t -fit) 

This added to the second term gives 

(B26) (1 + <5¥) u° (t - fit) = u Q 

or 

(B27) V(fit) R° = R (or R„) 

O ^ 

or 

H 6 t) 5 ° = 3 0 < or 5 C ) 

The superscript "o" denotes an elastic field entity, while the subscript 
”o" denotes an initiation level of the same entity (a material constant). 
Since the propagation occurs here at a certain steady-state value of R, 
the initiation and the critical levels can be regarded equal. 

Equation (B27) is identical with the Knauss-Dietraann equation given 
for cracks moving in linearly visco— elas tic solids [A7]. A somewhat 
different result was obtained by Kostrov and Nikitin [A 8 ], but their result 
can be shown to converge to (B27) if one assumes neglegible plastic effects. 
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c. Comparison of Slow Growth Theories Based op (a) Cherepanov 1 s, (b) 
Morozov^ and (c) Wnuk*s Criteria. 

(a) Constant Fracture Energy Hypothesis (Cherepanov), 

Three hypotheses are discussed here. The first one was proposed by 
Cherepanov [2], The basic physical assumption made by Cherepanov is that 
the total work done in separating two surfaces during an incremental growth 
is a material constant. This statement expressed in terms of elastic 
field entities and with the assumption that the Dugdale model applies, reads 

(Cl) 4 / Y 5u (x, Q( 2) ,2)dx * 2 5 52 

2 c 
c rf 



or 

(C3) 2Y[-^ + -^ ^ ] J u[x,Q(2) ,2}dx + 2Yu( tip) = ^7 

(compare Wnuk [9]). The above relation describes a slowly moving crack 
within the sub critical range of the applied load Q. Symbol u denotes the 
displacement perpendicular to the crack surface, 2 and a denote the half- 
length of the crack and the half-length of an extended (Dugdale) crack, 
respectively. Formula (Cl) is valid for an elastic-plastic solid which obeys 

the Treses. yield condition, and it does not account for the rate Sensitivity. 

An extension of the Cherepanov theory was proposed for the visco-elastic 
solids by Wnuk [A3]. The governing equation of motion for such a case is 
only slightly different from (Cl), and it reads 

(C4) 2Y { / § ] dx) f(4/i) - 3 C 

Here u° denotes the associated elastic solution for the same boundary value 
problem, Y is the normalized creep compliance function, A is the characteristic 
dimension of the Neuber domain and 2 the rate of crack growth. 

Let us briefly present the essential results pertinent to the small scale 
yielding range. The integrals involved in (C2) or (C4) are evaluated as 
follows (R/2 « 1): 


J 


2 


9u 

D2 



3u 

3x 


u(2> - u (cl) = u(*) =(^)R(£,Q(£) 


(C5) 


Iq/» = 


4Y B z q,QU) « 2 3R 
W 3 ' 3 R 


DQ 


Combining the above results in accordance with (C2) and (C4) and recalling 
that "S can be replaced by 2Yu^ or 2Y(4 Y/ttE)R^ we arrive at 


(C6) 


K + 3 K di - R * 


for an elas tic-plastic solid, and 
(C7) (R + | R-ff-ff} 'H.n/b = R t 

for a viscoelastic-plastic solid. Note that the rate-sensitivity and the 
time-dependent properties of a solid are reflected in the latter expression 
by the creep function Y. Of course this function will depend on the propagation 
rate 2. To reduce the above forms to just one equation 
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which would contain only the sought-for function R=R(£) , let us eliminate 

■ 

the rates dQ/d£ and d£/dt. This is done in two steps. First ¥(A /£) is 
approximated by the first two terms of its McLaurin expansion 

(C8) ¥ = ¥(0) + ¥(0)A/£ , ¥(0) = l t ¥ (0) = B 

and then d£/dt is replaced by Q(dQ/d£) - ^. Since the R-function depends 
explicitly on £ and Q(£), we have also 


<C9) 

hence 


dR 3R _9R dQ 
d£ = "5X + 3Q d£ 


(CIO) 

and 

(Cll) 


— ^ — / dR jyi\ 1 3 R \ 1 
d£ " C d£ " W W 


¥ 


BA dR 3R, .3R.-1 

o M£ “ IF V 


Therefore, the governing equation of motion can be written for both cases as 

(C12) {R + |r - |& } * a - If) (-!§: rh = R* 

where the parameter C = B/(J describes the time- depen dent response of the 
material. Obviously, when C = 0, i.e. when either the rate-sensitivity 
B is zero or the rate of loading Q becomes infinite, equation (C12) reduces 
to the equation (C6) which describes propagation in an elastic-plastic, rate 
insensitive material. 

Let us illustrate applications of eq. (C12) for a case of a central 

crack contained in a large plate (plane stress only is considered). For 

this configuration the plastic zone size R is the following function of 
4*\jl 

crack length andj^loading parameter Q (- tt<o/2Y, & =: applied stress, 

Y = yield point^) 
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(C13) R = ^JtQ 2 a) 

2 2 

Normalizing R and £ by R (=ttK /8Y ), and inserting 

* c 

(C14) |£ = AQ and = ^Q 2 


9Q h a a 

into eq. (C12) produces 
f ■ — /dp . 

3 p ' , dC C' 


(C 15 ) { p +‘4 p&- £)}*{ 1 + c (r ) (4~- -) (2pO 1/2 } =1 


R* 'dc 


(C16) 


da _ 2 ~ Q 2 C 

dC Q2 C (2 QC+C ) 


3 " o 

Here p=R/R^, £=&/R a and C o =(B/Q) (A/R^) . The functions to be determined 
from the above non-linear first order differential equations are 


R = R(£) 


from eq. (C15) 


and 


Q = Q(&) from eq, (C16) 

We do not have an analogous treatment for a large-scale yielding rate- 
sensitive problem. However, if the rate-sensitivity is absent, the governing 
equation of motion, i.e. eq. (C2) , can be shown to take on this form (central 
crack configuration) : 

(C17) ^ = 2 ^1 ~ g (Q tan Q + £n cos 0) 

^ S ec2 Q - tan Q] 


for a plane crack, and 


(C18) 


C 2 q 3 /3 • dq , , , 

(1 - q 2 )*' d£ + (l-q2)*/^ 


1-d-q 2 ) 1 ^ 2 


= 1 


q =tQ 


for a penny-shaped crack. Of course both the above equations reduce to a 
common form within the small scale yielding range, namely 


M , 3 2-Q 2 C 

d£ 2 q 3 c 2 


(C19) 



(b) Criterion of Constant COD at Instability (Morozov). 


In 1969 Morozov [A2] proposed that the subcritical growth ought to be 

controlled by the COD criterion rather than a constant specific fracture 

energy concept. Since in the subcritical range the tip displacement cannot 

be considered constant (in fact it is a monotonic function of the load 

remaining in equilibrium with a crack of a given length) , one would only 

require that the final value of the tip displacement, attained at the terminal 

instability point, should be a material constant. This is exactly what 

Morozov suggested. In order to ensure a constant COD at the end point of the 

stable growth, he modifies the fracture energy S' in such a way that the 

c 

terminal instability is always reached in accordance with the COD criterion. 

Thus Morozov’s becomes a function of geometry and the current crack 

length. Let us denote this function by , where the letter "R” stands 

R 

for "resistance". A similar concept was proposed in 1968 by Wnuk [A4] and [A5]. 
The function "S is evaluated from the following equation 





dx 


(C20) - 2Y + -J- f u(x,l) dx} 

* « * «crit 

The expression contained in the bracket becomes a function of l only, since 
Q is eliminated through the use of the COD criterion. (This latter criterion 
is applied to evaluate Q cr £ t first^. The next step, according to Morozov, is 


to apply the energy balance equation which reads simply 
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Note that his operator 6 implies an identical operation as the one indicated 
in eq . (C2) . 

To summarize briefly Morozov's theory let us write down the essential 
results. First, the function ^^(£)for a plane crack is found to be 

(C22) {[1 - exp (-2/e)] 1 / 2 eexp(l/e) cos 1 [exp(-l/e)J - 1} 


and for a penny-shaped crack (compare Wnuk [A4] and [5])it is 

(C2M 5 { > 

Note that both expressions reduce to just Irwin *s for a sufficiently 

c 

long crack (£>>1)* Combining the above formulae with Morozov’s criterion 
for crack extension (C21) gives 


(C24) 


dQ _ 2{F(Q~t(QtanQ + & n cos q)} 
dC C 2 [q sec2Q] 


for a 2D crack, and 

ia - 3(l-q 2 ) 3/2 { G(C) - C [l-(l-q 2 ) 1/2 ] (l-q2,- 1/2 } 

«" ** ’ q-f Q 

for a penny-shaped cracks The geometry dependent functions F(0 and G(c) 
are given by the expressions inside the brackets of (C22) and (C23) , respectively* 
No rate effects are included in either of the above equations. When the 
plasticity effects are neglegible ( C > > 1 or Q-K)) , one can shew that both the 
equations (C24) and (C25) reduce to a common form 


(C26) 


iQ. _ 2 - co 2 
c 2 q 


This equation, surprisingly, has a closed form solution which was net mentioned 
by Morozov, and it reads 
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(C27) Q( S) = C 1 [4(? - c 0 ) + 

Obviously , the Initial condition, Q(4 o ) = Q q is included in the above expression. 
Due to the simplicity of eq. (C27) we may directly compute the load and the 
crack length at the terminal instability point. Equating dQ/dC to zero we 
get 

<C28) q crlt . /2/t crlt 

^crit = ” "2 2 o Qo 

where R Q is the plastic zone size at the onset of slow propagation and R A 

is the value of R at the fully developed steady-state growth. Formula (C28) 
appears to have some experimental support, cf. Sullivan and Freed [A.6], 

We may add that somewhat unexpectedly the governing equation (C26) may also 
be derived directly from the classical Griffith energy balance equation 

(C29) <5U = 4y<5 l 

if the operator 5 is applied in the following sense 
(030) «D (ll o + m U ) St 

With U = <o^/E, Q=tt ^/2Y and B.* = ttK^/ 8Y^ , the energy balance equation 

c 

(C29) supplemented by (C30) reduces to eq. (C26) . 

To conclude the section on Morozov’s criterion of fracture, we add 
visco-elastic properties into his model. For this purpose consider the small 
scale yielding range only, and apply Wnuk’s [9] result for the "effective 11 
strain energy due to a crack contained in a visco-elastic solid 

(C31) U eff = U° F(A /l) 
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Here, the associated elastic strain energy U° is (l/2)R* "fj and the 

c 

■ * t 

function ¥ can be approximated by 1 + ¥ (0) A/A , while % is replaced by 
R*Q(dQ/dc) Then, the balance of energy criterion 


(C32) ¥ (A/A) 6U° 

becomes 


2 5 SI 
c 


(C33) 


(Q 2 £ + Q £ 2 ||) = 


This is a non-linear differential equation defining dQ/d£ in terms of Q 
and C> so that the numerical integration by Runge-Kutta method presents no 

•it 2 

problem. If we omit , however, the terms containing (dQ/d$) , then the 
above equation simplifies to a form remarkably similar to eq. (C26) , namely 



Here, the constant C Q encompasses the following group of parameters 


o " Q ' R* 

The ratio A/R* can be roughly approximated by the quotient R q /R a , where R q 
denotes the threshold value of the plastic zone size. Obviously, with 
C o =0, we recover the rate- insensitive equation (C26) . 


* The test run on IBM 360 shows no appreciable difference in the shape 
of the integral curve resulting from the complete equation (C33) and the 
simplified one (C34) . 
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The criterion of fracture proposed in this work and in [8] requires 
that the amount of deformation which occurs within the process zone during 
the time interval just prior to fracture is a material property. The 

deformation involved here is identified with the increment of displacement 

♦ 

AuLp produced at the control point P, within the time 6t = A /£. Thus we have 
(C35) A Up = u(Xp,t) - u(Xp,t - <St) = R q (4Y/ttE) 


where the right hand side contains the material constants only. If the 
viscoelastic displacements are substituted into the above expression and 
R is considered small vs. the crack size, the equation (C35) takes on the 
form (compare eq. (B23)): 


&)g.nf4R/A)+ RQR/H) 
A + RCA (3R/aQ)~l 


Here C = B/Q , and B ( ='i r (0)) denotes the rate sensitivity. In order to adjust 
the right hand side of (C35) to the boundary condition at R(£=£ cr ^ t ) =R*> 
we have made an additional assumption regarding R q , namely 


(C36ct) R q = j g.n(4R*/A) 


With this substitution and upon normalizing both R and Z by the s teady-state 

R* 

value of the plastic zone sizej^eq. (C35) becomes 

(C37) ^ = h n(l/p) + Cp(3p/'»C) (9p/9Q)~ 1 

dS 1 + Cp(3p/3Q)-1 

p = R/R A , c = Jt/R* 

To apply this equation for a certain crack configuration, one has to supplement 
it with the function R(Q,&), which is known from the stress analysis of the 
corresponding Dugdale crack problem. Let us illustrate this point for the 
simplest case of a central crack in an infinite plate under tension €T , 
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We have here 


(C38) 


R = •“ A(ir G-/2Y) 1 , or , 
P = \ CQ 2 , Q = 7r er / 2Y 


Hence we find the derivatives required in (C37) 
(C39) fljr- /2 pF , l^.-p/C 


and we have 


(C40) 


d£_ = k an (1/p) + Co 2 /cJ2£? 

dC 1 + C^/j2^ 


as the governing equation of the problem. It can be readily transformed into 
the (Q»e) plane. With the use of (C38) we obtain 

_ in (2/q2 0 - Q 2 


(C41) 


dQ 


QCC2 + CQ) 


This is integrated numerically for a chosen set of the initial crack length 
and the initiation (or threshold) load Q q . The curves generated in this 
way, for various rate-sensitivities C, are then compared with those which 
resulted from Cherepanov's and Morozov's hypotheses. The graphs are shown 
in Fig. A3. One may observe rather pronounced differences between the three 
theories described here, especially for the parameter C approaching zero 
(i.e. fast loading, or neglegible rate-sensitivity) . For large C all three 
curves seem to converge. In the limit of infinite C the graphs reduce to 
a horizontal line, which corresponds to a "creeping crack", described in 
more detail in [9]. 
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D. " Resis tance" Function or Apparent Surface Energy . 

A brittle solid fractures when the elastic energy release rate ^3 

attains the critical level, say ^ , which in turn is related to the material 

c 

+ 2 v 

toughness and the Young modulus ( j = K /E ,V The Griff ith-Irwin criterion 

c c y 

for brittle fracture simply requires that the equality holds * 

This simple view complicates a little when we try to describe fracture 
in inelastic solids such as elastic-plastic, or viscoelastic-plastic materials* 
There the criterion for fracture (equation of motion) reads * 


®i> R < 1 + - R * 


This equation can be cast into a form resembling closely the classical 

criterion S » if we agree to modify its right hand side. Note first 
c 

that the expression (R/R^) is identical with the elastic energy release 

rate . Now we can re-write (Dl) as follows 


CD2) -fn l+|||f ) T (4/i) - \ 

or, finally 
(D3) 5 = \ 

where the new function is defined by (D2) and (D3) 

R 

<D4) "S K = ^ 1 + ^f f 1 ,(A/J)rl 


i 


*) To focus the attention we consider here only the results based on 
Cherepanov f s hypothesis, therefore eq. (Dl) is identical with equation (C7) . 
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The function will be referred to as the "resistance" function, or as 
K 

the "apparent surface energy." For an ideally brittle solid ^ » but 

otherwise it is a geometry-dependent and a rate-dependent entity. We believe 
that an experiment designed to measure the fracture toughness will in fact 
record the "apparent toughness", related to by the well-known formula 

k r-< e ^r> 1/2 ' 


To illustrate possible applications of the proposed equation (D4) , 
let us consider a central crack contained in a large sheet subjected to 
a remotely applied tension & . Let us consider first an elastic-plastic 
solid. The length of the plastic zone is then R = (l/2}£ Q , where 
Q = 7r<S/2Y and the derivative dQ/d£ can be found from the eq. (C16) 


(D5) 


d£ 


-3 2- 


Q 2 c 


Q3 e 2 


-} R* 


-1 


C - t/E, 


(the rate sensitivity is zero for an elastic-plastic solid while f = 1) . 
Inserting (D5) into (D4) gives the resistance function 


(I>6) 




Q 2 (0 


t 


where the function Q = Q(0 has to be determined by integration of (D5) * 

The curves ^ vs* £ will depend on the initial crack length and the threshold 
R 

level of the loading parameter Q, as shown in .Fig. A4« A similar equation 
for a crack opened by a pair of wedging forces is discussed in [A3]. For 
this configuration the equation analogous to (D5) is 


(D7) 


dQ 

dl 



2 7 

3tt g f 2 tt \ 

2Q Q2 


1] 


Q = TIP/2YR* , R = (R a /2tt 2 )Q 2 /c 
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while an analogue of eq. (D6) is 
(D8) *5 = 5 Q 2 /2 ^ 2 c 

K C 

To conclude this section let us take a look at a visco-elastic solid 
described by a Voigt element (E2t^2^ connected in series with a spring 
E^. If all plasticity effects are neglected, equation (D4) gives 

(D9) " ^/V(&/1) 

which for our model will read 

(DIO) = 5 c (1 + [1 - exp (-A/£t 2 )] } 1 

This defines the rate dependence of the apparent t often measured in 

experiments which involve a crack propagating through a visco-elastic medium. 

The graph constructed according to (DIO) is shewn in Fig. A5a. 

The above relation can be generalized for any visco-elastic solid, if 

instead of vs. the propagation rate one plots Z) vs. the current crack 
K R 

length. Then* for a central crack configuration, it may be shown that 

• 2 

T(A/£)-( Q . ) (& /&) . The resistance function is proportional 

Griffith applied ° 

to the reciprocal of ¥(A/i) , therefore 

on) 5 r - 5 c <> 

where n is the square of the ratio of the applied stress to the short-time 
Griffith stress. A graph resulting from (Dll) is shown in Fig. A5b. 



Part II. Dynamic Crack Represented by the Dugdale Model 

by 
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SYMBOLS 

= Dilatational wave speed 
C2 = Shear wave speed 
k - c 2 /c x 

c = Rayleigh wave speed 

s 

k s " c s /c l 

p = mass density 

y = Shear Modulus 

v * Poissons Ratio 


Y = Yield stress of material 
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ATKI NSQN - BRO BERG - DUGD ALE MODEL 
STATEMENT OF PROBLEM 

A crack is expanding in an elastic solid under the in- 
fluence of an applied tensile stress P at infinity. As shown 
in Fig. 1 a thin region ahead of each crack tip is deformed 
plastically, i.e., the stress in this region is assumed to 
be equal to the yield stress Y of the material. The length 
of the crack at a given instant is 2a and the length of each 
plastic zone is w = c-a. The crack tips are moving at a con- 
stant velocity d and the ends of the plastic zones are moving 

with velocity c. p 

rtimi 


X. 


Fig. 1 

The solution may be obtained by super-posing the problem of a 
plate with no crack under uniform tensile stress P and a crack 
expanding under the influence of pressure P over the portion 
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of the crack surface defined by j x | <a and a pressure (P-Y) 
over the surface a<|x|<c. Since the problem possesses sym- 
metry with respect to the x-y plane an equivalent problem 
for the upper half-space which shall be denoted as Prob. A 
may be formulated (see Fig. 2). Retaining consistency with 
Atkinsons solution the formulation is for plane strain. The 
plane stress solution differs only by a constant and will be 
obtained from the plane strain solution. 


PROBLEM A 


A pressure P is acting on the infinite strip |x|<at on 
the surface z s 0 of a semi -i nf i ni te solid z>0. A pressure 
(P-Y) is acting on the infinite strips at<|x|<ct of the sur 
face z = 0. The vertical displacement u z at z = 0 is zero 
for | x | >c t . 

Find the normal stress at the surface z * 0 for |x(>ct 

and the displacement u for |x|<ct. 

£ 

. i 


k— c=dt — * 





Fig. 2 
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In order to solve Problem A, the following problem is first 
solved (see Fig. 3). 

PROBLEM B 


A pressure P is acting on the infinite strip j x | <a t of 
the surface z = 0. A pressure q acts on the s tri ps at< | x | <ct 
whereas no pressure acts on the surface z - 0 for |x|>ct. 
Find the displacement u z at the surface z = 0. 

* 



Fig. 3 

In Fig. 3, c : is the dilatational wave velocity. 
SOLUTION TO PROBLEM B 

The boundary conditions on z = 0 are 


° z = -P 0< | x | <at 

= -q at<|x|<ct 

* 0 | x j >ct 


( 1 ) 
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For the problem where a pressure p* is applied to the strip 
0< l x l£^ 1 t anc * no pressure is applied for Jxj^jt the boun- 
dary conditions are 


o z - -p* 0<|x| > <a I t 

= 0 | x | >a 


( 2 ) 


and the following expression has been obtained by Broberg 
for 3 2 u z /3t 2 on the surface z = 0 (see Brobergs e q . (16) where 
his symbols q, 3 , x are equivalent to p*, eij/Cp and-Cjt in 
this paper). 


& 


- jeL J 

2-TTk 4 C?lXl W 




(3) 


where 

F(e 2 ) 



(4) 


An expression for 3 2 u z /3t 2 that corresponds to the boun- 
dary conditions of eq.(l) may be obtained by superposing 
three problems: positive pressure P acting on a strip ex- 

panding with velocity a, positive pressure q acting on a 
strip expanding with velocity c, and negative pressure q 
acting on a strip expanding with velocity a. This is shown 
i n F i g . 4 . 
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SOLUTION TO PROBLEM A 

Now use the solution of problem B to set up an integral 
equation solution to problem A. First cons i der an incremental 
load on the half space defined by the following boundary con- 
ditions (see Fig. 5) 

dc? z = 0 | x | <ct 

= -q 1 ( v ) dv ct< | x | <vt (6 } 

= 0 vt< [ x | 

where, 


or 4 


( 7 ) 
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. i 

* or 

« — c*t — *i 

* 



n 

i 


n 

.is 



- ^'(r\r) $r\T 


Fig. 5 


Making use of eqs.(l) and (5), with appropriate change 
in variables, the solution for 3 2 u 2 /3t 2 corresponding to the 
boundary conditions of eq.(6) is: 


it' 


_ fytnr) Anr 




i«|f(S)] *(?-&) 


- TUS c In 

zt 



( 8 ) 


In order to satisfy the stress and displacement boundary 
conditions for Problem A, superpose stresses do z over the 
range of eq . (7) together with the stresses 


f-P 



0< | x | cat 
at< | x | <ct 
ct< j x j 


(9) 


where q = P-Y, on the surface z » 0 in such a manner that 
3 2 u z /3t 2 * 0 for | x | >ct. 
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The superposed expression for 9 2 u z /3t 2 is obtained by 
adding (3) with q=q Q and the integral over the range c<y<c l 
of eq . (8 ) . Then 



Setting 3 2 u z /3t 2 = 0 for |x|>_ct and making the following 
transf ormati on of variables, 

— 7 - 7 = r * ^7 = s » 7-q'(v) = *(s) , (11) 

c/t cj 

we obtain, for gt<|x|<c , or equivalently for 
(c/Cj ) 2 <r<l , 



+ — JJ F(>J , < V- < I 


( 12 ) 



or rearranging terms, 
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-IT Im. [H [ y _ 

Re [Fl'/r-jp 


: l)V M 


<tl£) 

I- 0 

^ - I'l 

(r- s)‘ 

X. 

(s-r)/r 



- 2J*j&e 
c 7 


2(p.a„) (d r ~ if) 

1 ^ C ' (- £) 


or 


_ InO I ( Vril 


Re 


Sfen< r - ?.)*> - /<*#- 




.i. f(i=JVO 0( S )6s = “fr * 2j P ^£(^) 
li I rr ire, n (»•- 

V 


(13) 


.-1 


< r< 1, 

'"l 

Eq - (13) is a singular integral equation which may be written 


C^Y-}^') +- 7 = <^ r > ) % < *"< l 

l/» ? ^ 


(14 » 


where 


aW = £) , 

Re[F(K-r C ' ' 


0 


, < v- < k l 

1 C, v 


4 k V(l-^- r ) (r , |M 4 * <r<| 

(r- 2KM v c ' ' , , 

(15) 


• Z. 


b(r> - r- c /,i 


(16) 




- 2 C 
Tf C, 


1 ' <fi(s)&s +• $ .?),. ^.(r £i2^{i7) 

TT J~jT TTC, 

\ f 


cs 
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The solution of the integral equation (14) may be written 
in the form (see page 331, Muskhel i shvl 11 i . Singular Integral 
Equations ) 


i. c6(v) - i - fil 


^ L. <2-xp(-(o(rQ 

(**- Q/ + bl fc-> 



where L is an arbitrary constant and 

1 


id (r) 


r 

1 • v* 


JLili 


(c 2 /^) 


s-r 


ds 


(19) 


<B> (r) 



( 20 ) 


o)(r) is calculated in Appendix I of Broberg and thus, 
exp[-w(r)] is written 




- 


o , 

(r- (r- ZV.'-y - 4k 3 /0'v)(k l ~v-) c ‘ 


(l - v-) r(r-Kf) 


(r- <% x ) J(r-Z\C-)* +lfck 6 (l-0(v--k t -) 


> k’x h ^ 1 


( 21 ) 



41 


where k s is the ratio of the Rayleigh wave speed c s to the 
dilational wave speed c . 


Recalling from eq . (18) that g ( r ) is dependent on 0(r ) , 
it is seen that eq . (18) does not yet solve explicitly for 
0(r). Thus, the next step is to consider the integral on 
the right hand side of eq . (18). Denote this integral by J. 


J = 


( s- 9^,0 EtL>(s)3 bis) <%,(s) ^ 

(s - O / aX*'H-b 2 (s} 




( 22 ) 


and rewrite g(r) in the form 


g(r) = B + A/O-taVc, 1 )] 


(23) 


where 


_ 2 ( P-q n )a (a 2 - c 2 ) 


TT C 


(24) 


3 - 

TT C, 




TT C.. 


c x 


(25) 


Then J may be rewritten as 




or 




( $ - &x^p[(o(s )3 b(s) [ ~ +* ^ , 

(s-r) J CL 1 ( s ) 


J = g (r ) I j (r ) + [A/((a 2 /c \ )-r)] I x (a z /c J ) 


(26) 
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where 


i,(fi 


'c 

H 5 - 


( s ~ C^( S )1 b(s) 

(s-r)J a!i7) v lefts} 


& 


(27) 


since 


q(s)-Ol(\t) A r - — — 4=; — — ■ ■■ ft - 

s'-r = —»-[(*- <^ c .) (r-^)J ( s .a^)(y..*£J 


I t (r) and I A (a l /c^) may be calculated by contour inte- 
gration as in Appendix 2 of Broberg the result being:- 


1M = 


TT ( M ” *0 


^ < v < k* 


tt(h A , n4kVir-k, t )(i-v , )iv-k^(r-^;) ^'<,<1 ( 28 ) 
(r-Zk^t )Gk 6 (l-v-’)(>'-k z -) 


and 

I^iVcJ) = tt[M- (a 2 /c 2 )] - Tr[(a 2 /c 2 )-k 2 )T (29) 

where 

( “£ f - Ik'*)’' - 4k7(>- <?% r Xk*- ^0 (30) 

with M being a constant it is not necessary to evaluate. 

Upon substitution of the values for I x ( r ) , I (a 2 /c 2 ) 
into J(r) and then into eq.(18) the resulting expression for 
<f>(r) is: 
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v(r-k^) 


A(m-v) ^ A ( M- 

v-(*--W, l )(v- $j») v "('" k s'l ( r - a /cf) 


Hr-V?) 


A(M -*/x) 1, 

v(>-- 


(r-Zk^- AkV(wvXk-r) 

Ji^iFW 



<V-<k l 


tr- 2 ^ , k l <v< I 

yci-vK >--%>? 


(31) 


Next* isolate coefficients of the term l/(r-k*) in eq.(31), 
thus rewriting <J>(r) as: 




-BM4-A + TA +L 
W/Cv'-W/') 


^BM-A-TA-L 
V" k* 


( 





( V-Zipf- 4kV(l-v-Xk x -y-' ) 



zrzk 1 - 




, k l < v-< | 


(32 ) 


The first term on the right hand side of eq.(32) does not 
satisfy the displacement boundary conditions for the problem 
(see Broberg). Then in order that this term be zero, set 


L = BM - Bk* 


A - TA 


(33) 
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and, is then, 


0M = |*B- 


T A 


h &7 


- 4kV(l-'~Xk t -y) 

S'.) 3 

l *•/ Ci-ocv- ^ 




V 1 < v-< 1 


(34) 


In order to evaluate the unknown constant 3 substitute 
eq.(34) into eq,(25). The result of this substitution is 


BZ 1 =(2q 0 c/ Cl ) - [TA/(a 2 /cJ)]Z, 


(35) 


where , 

2. 


an d 


2 C ( tf-crt - 4 k l ( <&t) + +k 4 ] C, Kfr-^r N 
c(.-c^) K{f^) 

_ 2[ ( c 4y-4kYl^k t )( c 4' t )t8k 4 ]c, i t r / j "y g-,\ 


2, 


2 ( fc^K ^jjk 0 - 8k 4 ( % - fjcj) 

(c^)(c4 r ^)( t _c> c ) 


(36) 


E(F^) 


\; 
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-h 


(< 


2(&c,)( jj -/ (i-’kf) i, ^ 

i-^)(v^) '<-*y s 


_' gk a c- 


jrUiz^l 




(37) 


where K and E are complete elliptic integrals of the first 
and second kind and where'll is a complete elliptic integral 
of +he third kind. ( see Byrd and Friedman, Handbook of Elliptic 

Integrals). 

STRESS ON SURFACE Z = 0 


The expression for stress over the interval ct<|x|<vt 

on the surface 2=0 may be found by integrating eq,(6) over 

+V\£ 

the range of eq.(7) and recall i ng^transformati on of variables 
in eq . (1 1 ) . Thus , 




(*)=- 


(<vr) A«\T 



As 


(38) 


which upon substitution of eq.(34) reads. 



\ 



( As 
sys(i-s)(s-^.) 


3 





Jk l -$ As 
S(s-a>yy 5 ( s _c^ r 

" (39) 




J 


C't C\K\< C. 
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The objective next is to make the stress singularity at 
xj = ct equal to zero. This is accomplished if 


B = AT/[(c/ Ci ) 2 - (a/cj 2 ] 

with the stress expression becoming, 

cr M = - f 

. 2(c' l -o :')[ J s(s-^l) 


(40) 


tfi 2 - 


- 4-k ? f ^ 1 , t't<lx\< 

k*-**V*c.-<*> j (4i) 

where it can be seen by inspection of the integrands in eq .(41) 
that the integrals remain bounded when |x| = ct. 


DISPLACEMENTS 

The region of non-zero displacement on the surface z = 0 
is | x | < c t . For this range, eq.(10) reduces to. 


_ 

at 1 - 



+ 


(P-foV 


I _ ait 1 

x v 



nr c, 

I 

C" x 


ce. 


fy'(cv-) AflT ,2 = 0, ocjxKc't 


^ j 


( 42 ) 
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The change of variables in equation (11) will again be 
used. Thus, eq.(42) is rewritten in the form, 


d^Vj _ -*ReLF(^r)3 j feckyr (P- qj cU t r 

dt 1 ZiTk 4 ^c l L lxl^ (|r- c^) (>-_ ar> c ^) 


ML k-s 

v 

/ c. t 


7#^- .1 *(.) 4* 
\/r(ir - <Vc*)J 


') ^^WoJUavU^ ; '0 Or ( /"^V 


(43) 


Eq.(43 may be rewritten as 


A> s _ -»-'Re[KV 0 f _$_£! , ( v- Q.>) (P-O ol 

2nk^ixi(v- «|j.)L Ci ( *•- <Sh*l c, 


-if ^)As +±f L0(t)k 


0(r( (44) 


Further simplification is allowed through eq.(25) so that the 
preceding expression becomes 


3*w. _ -fTfeLFWI ttB (^-c'KP-gJa: 

■>+* Znl^e'*' 2- c 


^ 5 > As 




( 45 ) 
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It remains to evaluate the integral on the right hand 
side of eq. (45). Making use of eq.{34) 



The four integrals in the preceding expression will be evalu- 
ated by contour integration (see for example p. 276 of Muskheli- 
villi)* Note that all the following are evaluated for the 
range 0<r<(c 2 /c 2 ). 


i.- 


f ClzA k * - ) 1 

' C Z 5 ( S - V-yO-s)(s-C^) 


(47) 


The real integral I may be 
written as the limit of a 
contour integral in the t-plane. 

2.T - dlf 

l ' ' Jp r(+--rVCi--*>c t- 

= - 2 Tr(f T M- G-^Cv) - , °CV-< 

where 



f A (r) = - 


(r-2k 2 ) 2 

r/(l -r ) [ (c 2 / c 2 ) -r] 


Muskhel i shvi 1 1 i , N. I., Some Basic Prop! ems of tne Mathemati cal 

T h eory~~ of El as t i c i ty 
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I 4k*c 

<\< r > * - -Sr" 1 


G 1 ( r ) = 1 


>4 


Then , 


T v = tt( 


(v~-3k>V 


_ 4k z c. 


v/t'-v-XS^-lr} a*v- 


+ J 


) 


( 48 ) 


I,- 


As 


,. 1 V 


'c? 


1 


( k 2 - 1 ) d t 


21 „ = Tim 4 , 

r t(t-r)/('t-kO(t=c 2 /c|) 
2 

= -2ir|f n (r) - G J 1 ( r ) - G^(r)j 



where , 


f H (r) 


r— k ; 


r/lYa-r) (c*/cf - r) 


G”(r) = 


kc . 
3r” 


6l r (r) = 0 


\; 


Then , 
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= -Ztt( 
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= o 
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Using the preceding expressions for I, , I z , I 3 , and 1^ 
eq.(46) reduces to, 


f 0(s] As = TTB / (r-lW X -A\cy (t-<rW.*i 

J (s-V-'l V rjC\-' r) 0$ c> 

n ^ 


H* 1 




-uT/v{ 


4kVn-vKv^-^ 


( - 2 kf - 4 kV( I - ( |f - 4>,4 1 

( *■ ■ &V) /('" "nr ^ 


and 


a\j* - 


2tt£ k/Mxl ( v*") 


i r - 


_ HB (r-Zk^- 4kVl-rXk^r) 

( v- zgf - 4 kV(\ -V-) ( 
z ^ <*f - O 

_ . ~ jjgj - 4k^Xi^v)(k. t --^ t Q 


(56) 


| 0<r< c% fi7 , 

(«>c>)(r- ^/C'-^X^r Cf 


But, recalling eq.(4), 


■MfW> 


_ 4k 4 




(r-^k L l l - 4kJ(k l -vll-v) ' 0< '" < > ( 58 > 


and making use of eqs.(24) and (30), e q . (57) reduces to 
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/u, _ ! 3^- - X* »* 

£Ut c^-v-11^7 ( 


^ >•*• \ 

- J 


For the first integration note the following 

2 > Br 5 __ 2 -^ S - _ 7 C,r Vt - a 

5t Jr “ |x| 2v* 


r< ^ 

s 

( 59) 


Hence 


i thi)- — -f JL_ 

in 'it I 2^c,l r) 3/l * 


and 


3>u 


1 


3 


TA 


An 


3 + WJ( ^ SUJ (r- 45 jyspr 


TA 


C&V* 




Making use of eq. 212.9a of Integratafel , 


3l> £ _ 

ar 


2eq 
TA 


/c ' 1 a ' 1 J / c ' 1 V 


Hmf 


I /t' 4 _ r 

I vLc£ r V c. 1 c, , *~ 


(4'TT , AT i 


} 


(60) 


(61) 


(62) 


Recalling the finiteness condition as expressed by eq.(40), 
eq . (62 ) reduces to 


£_v % — 


= ■ -~ I£l 


Zee, (<%!.- a ?3, v ) /i ' 


jinn 


TIE 

Mr* *Mr% 


, OtK‘M 


( 63 ) 
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V 

Again writing the derivative in terms of r r a second 
integration with respect to r leads to the result, 


u - -^ TA 

* / .7 


fjL fa 

Mr 0 

M -Mcr r 

v-ku 

M r l 

+ JM~ r 




%■ + 


- t/Hr* 


OC \T< 

(64) 


In terms of x and t, the result of eq.(64) may be 
written 




tMZ - M\ x - ^ 


t x Qm 


x/cVc? - 


1X1 < c f 


(65) 


Half of the cracks opening displacement, 6y, is found by 
letting [x| -*• at. The result is 





( 66 ) 


It should be noted that the expressions for displacement 
make use of the finiteness condition, eq.(40). Then, by eqs. 
(35), (36), (37) and (40) the following relationship between 
a, c, P and Y is obtained. 
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(fc / 4 ) 



Fig. la . R-curves for a rate sensitive Tresca solid at various levels of the rate sensitivity /speed 
of loading ratio. Ultimate instability is reached when the slope of the R-curve plotted 
on log-log scale equals unity (marked by circles). Crack of initial length equal lOOx 
characteristic dimension of Neuber's domain is embedded in an infinite plate. 





iOfXdln^ 


■ lM 
IIM 

■ MB 

■ MM 

■ M ■ 

■ MB 

■ MB 

■ MB 


Fig. lc Slow growth at various loading rates. Loading parameter is 
plotted vs. current crack length at a chosen loading rate 
and material rheological sensitivity. Initial crack size is lOK^. 
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Fig* Id slow growth at various loading rates. Current crack 
length is plotted vs. time at a chosen loading rate 
and material rheological sensitivity. Initial crack 
size is 10 R . 
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Fig. 2. Sub critical growth of a crack embedded in a viscoelastic solid. 
Concentrated loads are applied directly to the crack surface. 
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Amount of slow growth Ac - Cr - C q vs. initial crack size at 
various rates of loadings and a constant initial crack 
length/panel width ratio. 
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Fig. 6. Amount of slow growth Ac = £f - £ vs. 

initial crack size /panel width ratio at 
a constant initial crack size and a 
variable rate sensitivity /rate of 
loading ratio. 
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Amount of slow growth A£ = t f ~ C 0 vs. rate sensitivity /rate of 
loading ratio. Initial crack size is 10,100 and 1000 in 
units in Fig. 7a. In Fig. 7b the Initial crack size is constant 
and equal 10R*. while the ratio p - (ir/2)*crack length/panel width 
varies. 









Applied load as a function of the increment of crack length 
generated during the slow growth stage which originated at 
various initial flaw sizes. Kate sensitivity C and initial 
crack size/panel width ratio are constant. 








Fig .A2. Front of an advancing crack and the associated yielded 
(or crazed) zone. P is the control point. 








ED APPAR 











CRACK LENGTH /INITIAL CRACK LENGTH, x =-t/4 0 — * 

FIG. A5 APPARENT SURFACE ENERGY (RESISTANCE) OF A 

VISCOELASTIC SOLID IN THE SUBCRITICAL RANGE, 
AND UNDER SUSTAINED LOAD, m = APPLIED LOAD/ 
GRIFFITH LOAD 




